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Abstract 

In quantum electrodynamics, static electric fields are screened at non-zero 
temperatures by charges in the plasma. The inverse screening length, or De- 
bye mass, may be analyzed in perturbation theory and is of order eT at rela- 
tivistic temperatures. An analogous situation occurs when non-Abelian gauge 
theories are studied perturbatively, but the perturbative analysis breaks down 
when corrections of order e^T are considered. At this order, the Debye mass 
depends on the non-perturbative physics of confinement, and a perturbative 
"definition" of the Debye mass as the pole of a gluon propagator does not even 
make sense. In this work, we show how the Debye mass can be defined non- 
perturbatively in a manifestly gauge invariant manner (in vector-like gauge 
theories with zero chemical potential). In addition, we show how the 0{e^T) 
correction could be determined by a fairly simple, three-dimensional, numer- 
ical lattice calculation of the perimeter-law behavior of large, adjoint-charge 
Wilson loops. 
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I. INTRODUCTION 



Electrically charged particles in a hot plasma react to electromagnetic fields and cause 
screening of static electric fields at large distances. The inverse screening length, also known 
as the Debye mass ma, may be computed in QED by considering the exchange of a single 
virtual photon between two static test charg depicted in Fig. [I|. 




FIG. 1. A single virtual photon exchanged between two static test charges. 

The long distance fall-off of the static potential is determined by the position of the pole in 
the photon propagator at zero frequency. This is given by the solution = —m^ to 

/ + noo(o,p) = o, (1.1) 



where n^i,(po,p) is the self-energy of the photon. In the ultrarelativistic limit (when particle 
masses and chemical potentials are negligible), the leading-order result is easily computed 
from the one-loop graph of Fig. |^ and yields 




FIG. 2. One-loop self-energy of a photon. 



For reviews, see refs. [|]J^. 
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FIG. 3. One-loop self-energy of a gluon. 

Unlike electric fields, magnetic fields are unscreened, which is reflected by the fact that 

lim n,,(0,p) = 0, 2,j = 1,2,3. (1.3) 

A similar perturbative calculation carried out in non-Abelian gauge theories, using the 
one-loop graphs of Fig. |^, yields a lowest-order result of 

md = mo + 0{g^T) , (1.4a) 

where 

"^0 = + -^j 9T (l-4b) 

for SU(A^) gauge theory with rif Dirac fermions. 

For the sake of better understanding the nature and reliability of perturbation theory at 
finite temperature, there has long been an interest in computing the leading correction to 
this result 0]. It is known, however, that this correction cannot be computed perturbatively 
in non-Abelian gauge theories As we shall review below, the 0{g'^T) correction to the 
Debye mass receives contributions from fundamentally non-perturbative physics associated 
with the interactions, at high temperature, of magnetic gluons with momenta of order g'^T. 
The best that can be done perturbatively is the extraction of a logarithm at that order : 

= mo + ^Ng^Tln (^^j + cg^T + 0{g'T) . (1.5) 

The constant c, however, is not computable by perturbation theory. 

Because the physics of the 0{g^T) correction is non-perturbative, it behooves us to 
formulate a non-perturbative definition of what we mean by the Debye mass in the first 
place. Such a definition should be gauge invariant and preferably implementable in numerical 



lattice simulations. The definition (|1 . 1| ) is unfortunate in both these respects. In particular, 
the self-energy is not itself gauge invariant in non-Abelian theories. There are formal 
proofs that the pole position is gauge invariant order by order in perturbation theory [Q, 
but this is of limited use since perturbation theory breaks down beyond leading order. We 
should look instead for a definition that is manifestly gauge invariant. 

The purpose of this paper is two-fold: (i) to give a natural non-perturbative definition 
of the Debye mass, and (ii) to show how the constant c in the expansion (|1.5| ) could be 
extracted from a relatively simple numerical computation of the perimeter law fall-off of 
large, adjoint-charge Wilson loops in three-dimensional, zero-temperature, pure lattice gauge 
theory. In section ||, we briefly review the source of the breakdown of perturbation theory. 



In section |T|, we construct a manifestly gauge invariant, non-perturbative definition of the 
Debye mass. We review why one method sometimes suggested in the literature — extracting 
the Debye mass from the long-distance correlation of Wilson lines — is inadequate. Our 
definition works only for vector-coupled gauge theories, such as QCD or QED, and only at 
zero chemical potential. We explain what the difficulties are for axially-coupled theories or 
non-zero chemical potentials, and we outline the problems with making a non-perturbative 
definition of the Debye mass in those cases. Finally, section |V| contains our derivation of 
the 0{g'^T) correction to the Debye mass in terms of three-dimensional Wilson loops. 

II. NON-PERTURBATIVE HIGH TEMPERATURE PHYSICS 

The physical picture behind our definition will be clearer if we first review the source 
of non-perturbative effects in hot non-Abelian gauge theories.@ The problem is easiest to 
understand by considering a series of effective theories corresponding to larger and larger 
distance scales in the hot plasma. Since we are interested in studying the screening of 
static electric fields, we can work directly in Euclidean space where non-zero temperature 



^More details can be found in the discussions and reviews of refs. |i|,^,^0|jTT| . 
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corresponds to making the Euclidean time direction periodic with period (3 = 1/T. So 

Z = J[VipViljVA]exp(^-^ J^^ drd^x C^^ , (2.1) 

where we have suppressed details of ghosts and gauge-fixing. Boson (fermion) fields have 
(anti-)periodic boundary conditions in Euclidean time. At distances large compared to (3, 
the dynamics of the time direction decouples, and one obtains an effective t/iree-dimensional 
theory of the zero-frequency modes of the original fields. Since the fermionic fields are 
anti-periodic, only the bosonic degrees of freedom are relevant in this effective theory.li 
Schematically, 

J [DA] exp (^-^ J d^x C,s^ . (2.2) 

The effective theory ( p.2[ ) is a three-dimensional gauge field A coupled to a three-dimensional 
adjoint-charge scalar corresponding to Aq^ One of the effects of integrating out physics with 
momenta of order T is that the adjoint scalar obtains a mass of order gT. As indicated in 
( ^^ , the gauge coupling constant in the three-dimensional theory is g"^ = g^T. 

Next, consider distances in the three-dimensional theory that are large compared to 
1/gT. At these distances the adjoint scalar decouples, and the new effective theory is a 
pure gauge theory in three dimensions.! Three dimensional non-Abelian gauge theories are 
confining. Moreover, the only remaining parameter of the theory is the three-dimensional 
coupling g'^T and so, by dimensional analysis, the confinement radius is of order 1/g'^T. The 



•^Physically, light bosons dominate over light fermions at low frequency because of the infrared divergence 
of the Bose distribution l/(e^-^— 1) as £'^0. 

^Purists may object to saying that the adjoint scalar in the effective three-dimensional theory arises from 
Aq in the four-dimensional theory, because this statement is gauge dependent. The gauge independent 
identification is that the three-dimensional scalar corresponds to the traceless part of the path ordered 
exponential Vexpi J^^ Ao{t, x) dr in the original four-dimensional theory. However, we shall continue to 
refer to the adjoint scalar simply as Aq. 

^There can be scalars too if they are part of the fundamental theory, such as for electroweak theory, and 
if the temperature is fine-tuned to be near a phase transition. 
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physics of magnetic gluons with momenta of order g^T is therefore non-perturbative, and the 
physical states of the three-dimensional effective theory are glue-balls rather than individual 
gluons. This is unlike the case in zero-temperature /our-dimensional theories, where the 
confinement radius diverges exponentially as (7 and non-perturbative contributions are 
never the same order as perturbative ones. 

It is important to keep in mind that the physics at large distances is the physics of three- 
dimensional confinement, and that it is this confinement which cuts off infrared divergences 
encountered in perturbation theory. Some papers in the literature work under the misappre- 
hension that the infrared physics is instead cut off by some sort of mass of order g^T for the 
gauge field A. This is as misleading as thinking of confinement in zero-temperature QCD as 
being described by a gluon mass. A mass would cause large (fundamental-charge) spatial 
Wilson loops in high-temperature gauge theory to have perimeter-law behavior because it 
would screen the gauge force. Instead, such loops will have area-law behavior — the signal of 
confinement. 

III. DEFINING THE DEBYE MASS 

Begin by considering QED. The Debye mass can simply be defined by the correlation 
length of the equal-time electric field correlation function: 

{E{x) ■ E(0)) ~ e-^'"<^l3;l/|^|3 as \x\ 00, (3.1a) 

or 

md = - lim \x\'Hn{E{x) ■ E{0)) . (3.1b) 

This is equivalent to the definition ( |1 . 1| ) in terms of the photon self-energy because the 
exponential rate of decay of a propagator at large distance is determined by the location of 
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singularities nearest the real axis in momentum space.Q 

Unfortunately, this is a poor definition in non-Abelian theories because E is no longer 
gauge invariant. One might instead consider a definition in terms of the correlation between 
two static test charges. Specifically, a manifestly gauge-invariant possibility would be to 
define the Debye screening length as the correlation length between Wilson lines (also known 
as Polyakov loops), 

{L{x)L\0)) I e-"^l^l/|s| as |a;| ^ oo , (3.2) 

where the Wilson line 

L{x) = tr V exp i / AQ{T,x)dT (3.3) 
Jo 

is the trace (in the fundamental representation) of the path-ordered exponential of the line 
integral of the gauge field around the periodic Euclidean space. (V denotes path ordering.) 

Although this definition has occasionally been suggested in the literature, it is wrong. 
Even in QED, it fails to isolate the quantity one wants to identify as the Debye mass. Though 
Wilson lines couple directly only to electric fields, they couple indirectly to magnetic fields 
through interactions, and magnetic fields are not screened. Fig. ^ shows how two Wilson 
lines can exchange a pair of magnetic photons in QED, and so, despite the screening of 
electric fields, the correlation (0) falls off algebraically instead of exponentially.! In non- 
Abelian gauge theory, the coupling to the spatial gauge field can be even more direct [^0, 
as in Fig. |. The non-Abelian case is slightly different from QED, however, because three- 
dimensional confinement implies that the Wilson lines cannot exchange a massless pair of 



^The fact that ( |3.lD specifies coincident times while (f.l) refers to zero frequency makes no difference. 
All contributions from the (discrete) non-zero frequencies to the equal time correlation function decay as 
0(e~^'^^l~l) or faster; hence the zero frequency component dominates at large distance. 

^The potential between arbitrarily heavy test charges, V{x) — — /3~^ln (L(a;)Lt(0)), decreases as Isl"^) 
reflecting a magnetic Van der Waals interaction between the two electron-positron clouds screening the test 
charges. An analogous case of algebraic screening in non-relativistic theories at finite density is discussed in 
ref. il. 
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magnetic gluons; the pair will instead form a glueball with a mass of order g^T, and so the 
correlation length of Wilson lines defined by ( p.2D will be of order l/g^T Regardless, 
this is not the physics of electric screening. 



FIG. 4. Power-law interaction between two Wilson lines, representing static test charges, due to the 
exchange of a pair of unscreened, magnetic photons. 



Fortunately, there is a simple symmetry which can be used to exclude the unwanted 
exchange of a pair of magnetic photons or a magnetic glueball: Euclidean time reflection. 
Euclidean time reflection corresponds to what in real-time is called TC, or time reversal 
times charge conjugation,! and the crucial property is that Aq is intrinsically odd under 
this symmetry while the spatial gauge field A is even. The Euclidean description is more 
convenient for our purposes, and we shall frequently refer to the symmetry simply as "time 
reflection." (The reader should note that in Euclidean functional integrals time reflection is 
no more subtle a symmetry than spatial reflection; there is no extra complication associated 
with anti-unitarity.) Euclidean time reflection is a useful symmetry because, in the effective 
three-dimensional theory, the only effect it has is to negate the adjoint scalar Aq. 



mnemonic for this fact: In real time, CVT must be a symmetry of any Lorentz invariant (and unitary) 
theory. In Euclidean space, VTZt, where TZ^- denotes Euclidean time reflection, is a pure rotation and must be 
a symmetry of any Euclidean invariant theory. So TZr must correspond to TC, since V is time independent. 





FIG. 5. An exchange of magnetic gluons between two Wilson lines. 
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If one considers the correlation of a pair of time-reflection odd operators, instead of the 
Wilson lines, then the zero-frequency magnetic contributions of the type depicted in figs. ^ 
and 1^ will be eliminated. The lightest intermediate states which can contribute will be 
those containing a single Aq (plus surrounding glue), so that the correlation length will be 
mo + 0{g^T). Any local, gauge- invariant, time-reflection odd operator can be considered as 
a replacement for the Wilson line, leading to a general definition of the Debye mass: 

Definition. Consider the correlation lengths defined by the fall-off, at large spa- 
tial separation, of the correlation {A{x)B{0)) between operators A and B that 
are local (in 3-space), gauge invariant, and odd under Euclidean time reflection 
{i.e., real-time TC). The inverse Debye screening mass l/m^ is the largest such 
correlation length. 

We are thus able to define the Debye mass directly in terms of the long-distance fall-off of 
certain correlation functions. This definition will only work, however, in theories where real- 
time TC is a good symmetry; otherwise, there is nothing to prevent states with a single Aq 
from mixing with A glueballs, and all of our inverse correlation lengths will again be 0{g'^T) 
instead of 0{gT) and will be unrelated to the physics of electric screening. The restriction to 
TC-conserving theories means that the Debye mass cannot be rigorously defined by the long- 
distance fall-off of correlation functions in theories with axial couplings, such as electroweak 
theory, or in the presence of a non-zero chemical potential. We shall comment again on 
these cases later, but for now our discussion will be restricted to vector-coupled theories, 
such as QCD, at zero chemical potential. 

Before proceeding further, it will be convenient to rephrase our definition in alternative 
language. Suppose that the separation x of our operators is in the z direction. In Euclidean 
space, there is nothing that distinguishes the time dimension as fundamentally different from 
the spatial ones. One may turn one's head on the side and interchange the labels z and t, 
as depicted in Fig. |^. One then interprets the original four-dimensional field theory as a 
zero-temperature theory with one periodic spatial dimension, instead of a finite-temperature 
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FIG. 6. Interpreting Euclidean time as a periodic spatial direction, by relabeling coordinates (r, z) as 
(z-,f). 

field theory with all spatial dimensions infinite. Our correlation functions are now correlation 
functions with large separations in "time," and their exponential fall-off is determined by 
the energies of the physical states in this zero-temperature, spatially-periodic field theory. 
So the following is an exactly equivalent definition: 

Alternative definition. Recast the theory as a 3+1 dimensional field theory at 
zero temperature, where one of the spatial dimensions — call it z — is periodic 
with period /3. Then, in a Hilbert space interpretation, the Debye mass is the 
energy of the lightest state that is odd under z refiection. 

For future reference, we should clarify that we will always use the notation Aq to denote 
the component of the gauge field in the periodic direction, regardless of whether we are 
interpreting Euclidean "time" according to the original definition or the alternative. 

In addition to their behavior under z reflection, henceforth denoted TZ^, eigenstatcs of the 
spatially-periodic 3-1-1 dimensional theory may also be classified according to their behavior 
under other space-time symmetries. Specifically, excitations at rest {px — Py — 0) can be 
assigned quantum numbers J*", \pz\, and Rg where J is the angular momentum in the xy- 



T 
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jC(P) 


Rz odd operators 


0-+ 


tr (Imfi) = Im L 


or 


tr(Fo3 (Fi2)2) 


0+- 


tr (ImfiFia) 


or 


tr (Fo3 F12) 


0— 


tr(Imfi[Fi2, (F,3)2]) 


or 


tr(Fo3 [F12, (F,3)^]) 


0++ 


tr(Imfi[Fi2, [Fi3,F23]]) 


or 


tr(-^03 [-^12, [-^13,-^23]]) 


1+ 


tr {ImQ F,3) 


or 


tr (Fo3F,3) 


1- 


tr (Imfi{Fi2,F,3}) 


or 


tr (Fo3 {F12, Fis}) 



TABLE I. Examples of gauge- invariant Euclidean time- reflection (or TZg) odd operators which couple to 
specific J'-^^-^^ sectors. Here, i = 1, 2 is a 2+f dimensional spatial index and Q. — V exp dx^ denotes the 

un-traced Wilson line or Polyakov loop. ImTl is shorthand for the anti-Hermitian part, ImSl = (O — ^^^)/2. 

plane, C is charge conjugation, is the momentum in the periodic direction and R2 is the 
sign acquired under z reflection. For J = 0, there is also one additional quantum number: 
the sign P of the state under two-dimensional reflections.!^ The lowest energy states will 
have = 0, and our definition of the Debye mass restricts us to = — ; so the states of 
interest can be summarized by J^^^\ It is not clear a priori which J^^^^ sector will contain 
the lightest R2 odd state. 

Gauge invariant operators which couple to specific Rz odd symmetry channels may be 
easily constructed. For example, under TZ^ reflection Aq —>■ —Aq and L L'^ . Hence, the 
time-reflection odd part of the Wilson line is just the imaginary part. ImL is also odd 
under charge conjugation, but is even under x- or y-reflections, so the imaginary part of the 
Wilson line probes the O""*" sector. Table H illustrates some of the possible gauge invariant 
operators which can be used to probe various symmetry channels. In the language of the 
effective three-dimensional theory (now to be regarded as 2-1-1 dimensional), each of these 
operators creates an Aq accompanied, because of confinement, by a neutralizing cloud of 



^The representations of 0(2) — Z(2) x S0(2) are (i) two-dimensional representations for each non-zero 
value of J^, and {ii) two one-dimensional representations, distinguished by their Z(2) charge, for — 0. 
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glue. The lowest mass state in each of these channels will have a mass of tuq + 0{g^T). In 
a direct lattice determination of the Debye mass, one should in principle check these, and 
perhaps other, channels in order to find the lightest state. (Alternatively, one could consider 
correlations of operators with less symmetry.) 

In pure SU(2) gauge theory,0 the operation of charge conjugation is in fact an element of 
the gauge group (namely ia2)- Hence, in this theory, any gauge-invariant state must have C 
even, and so the possible sectors are restricted to the J"*"^^-* channels. (Note that the charge 
conjugation odd operators shown in table | vanish identically for SU(2), as they must.) 

Wilson lines in pure gauge theories 

In pure gauge SU(A^) theories (that is, gauge theories without matter field^ll) there 
is one additional subtlety which occurs with the Wilson line ImL, or with more compli- 
cated operators containing a Wilson line wrapping around the periodic r direction. A 
Euclidean pure gauge theory at non-zero temperature is invariant not only under periodic 
gauge transformations; it is also invariant under non-periodic gauge transformations that 
globally multiply the fundamental representation Wilson line L by an element of the center 
of the gauge group.lii For SU(A^), the center is 'Zi{N), the A^-th roots of unity. This Z(A^) 
symmetry is spontaneously broken at high temperature, and there are A^ different (pure 
phase) equilibrium states distinguished by the phase of the Wilson line, arg (L) = 27Tk/N, 
k = 0, . . . , N—1. Only one of the equilibrium states (the one in which (L) ^ 1) is invariant 
under the naive definition of time reflection. (Each of the other A^— 1 equilibrium states 



^°These comments also apply to the pz=0 sector of SU(2) theories coupled to fermions (but not complex 
scalars), since the fermions are irrelevant. In addition, they apply to any group which, like SU(2), has only 
real (or pseudo-real) representations. 

^^Or more generally, theories whose gauge group has a non-trivial center but all of whose fields transform 
trivially under the center. 

-'^^For a review of this symmetry and its role at high temperature, see ref. H. 
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is invariant under a re-defined time refiection which combines the original refiection with a 
non-trivial 'L{N) gauge transformation.) 

In order for Im L to probe the Debye mass, one must work in the single pure phase equi- 
librium state which is invariant under (the chosen definition of) time refiection. Otherwise, 
time refiection will fail to select the charge-screening excitations of interest. However, a 
gauge theory functional integral which is invariant under the Z(A^) center symmetry neces- 
sarily averages over all N spontaneously broken phases. Because of this, the Im L correlation 
length, computed with a Z(A^) invariant functional integral will be 0{g'^T), and have noth- 
ing to do with the real Debye mass. This can be seen directly from the fact that the (LL) 
and (L^L^) correlations vanish by Z(A'") symmetry, and so ImL and ReL will have identical 
0{g^T) correlation lengths. 

This difficulty does not reflect any inconsistency in our deflnition of the Debye mass, 
because ImL is not actually gauge invariant under the full gauge group of pure gauge 
theories (which includes Tj{N) center transformations). Hence, in pure gauge theories, it 
does not meet the requirements stated in the deflnition. 

Nevertheless, one may avoid this difficulty in pure gauge theories, and obtain an 0{gT) 
correlation length for ImL, in either of two ways: change the operator, or change the 
theory. Fixing the operator is easy: simply replace the fundamental representation trace in 
the definition of the Wilson line by the trace in some other complex representation F which 
is invariant under the center of the group. For example, in an SU(A^) theory the symmetric 
tensor product of fundamental representations is suitable {i.e., the 10 of SU(3)). Then 



is invariant under Z(iV) transformations and is completely unaffected by the spontaneous 
breaking of the 'L{N) center symmetry. Or, one may use a local operator (of the same 
symmetry) which does not involve a Wilson fine at all, such as tr (F03 {F12Y). 

Alternatively, one may restrict expectation values to include only the single equilibrium 
state with {L)^l by adding an infinitesimal source to the Lagrangian that biases the system 




(3.4) 
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toward the desired Z(A^) sector :Ei 

= J[VtpV^VA]exp {-^j^drd^xC^ + e j d^xLj , (3.5) 



and then send e to zero after the infinite volume hmit. Therefore, the IniL correlation length 
^ may given by 

= — lim lim lim In (ImL(a;)ImL(0))^ y . (3-6) 

jI-^O e^0+ V^oo ' 

Adding the source term explicitly breaks the Z(A^) center symmetry, thereby reducing the 
full gauge symmetry of the theory so that L is now fully gauge invariant. Of more practical 
use for numerical simulations, one can simply run simulations in a large enough volume that 
one finds no jumping between the different Z(A^) equilibrium states, as measured by (L). 
One would need to take data generated by a run in a single pure phase where L=(trL) is 
non-zero, and then measure the long-distance fall-off of 

(Im[L(x)Z^]Im[L(0)I^]) (3.7) 

in lieu of (ImL(x) ImL(O)). 

Axial theories and chemical potentials 

As mentioned earlier, our definition of the Debye mass does not work for theories in 
which real-time TC is not a symmetry. Hence, it cannot be applied to gauge theories with 
axial couplings, or in the presence of a non-zero chemical potential. In the language of 
our alternative definition, where the z direction is viewed as "time," the problem manifests 
as follows: the lightest state with a single Aq is no longer stable against decay into an A 
glueball. Nevertheless, there is still a singularity in the complex momentum plane associated 
with this Aq "resonance." The situation is depicted in Fig. |^, where we have considered a 



"'^^The 7i{N) center symmetry is explicitly broken in theories with fundamental representation matter fields; 
hence in such theories no explicit symmetry breaking perturbation need be added. 
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FIG. 7. The singularity structure, in the complex pz plane, of a correlation of Euclidean time-reflection 
odd operators in a theory where real-time TC (a) is, or (&) is not, a good symmetry, tog is 0(g^T) and 
stands for the lightest glueball mass in the theory when the z direction is regarded as time. For simplicity, 
we have suppressed all singularities associated with excited glueball or excited Debye states. 

generic correlation of Euclidean time-reflection odd operators and sketched some features of 
its singularity structure in the complex pz plane. We have assumed p± = for simplicity. The 
case for a TC conserving theory is shown in Fig. 0(a), where the location of the singularity 
closest to the real axis is our Debye mass. Introducing a small amount of TC violation will 
mix in A glueball states, changing the analytic structure to that of Fig. 0(&). The Debye 
singularity from Fig. ^(a) is still present but has moved slightly off the axis onto the second 
sheet. 

One can still imagine, in principle, defining a (complex) Debye mass based on the location 
of the pole. This is, in fact, what one does every day when talking about the mass of an 
unstable particle such as the Z-boson or the 7r°. There is an important difference, however, 
which is that the z direction is not really a time direction and can be considered one only 
by analytic continuation. In the real world, one can reach the Z resonance experimentally 
by making —p^ + close to — M|. In contrast, one cannot experimentally study the static 
(po = 0) properties of the plasma by taking p^ close to —Tn\. And because the introduction 
of finite temperature breaks Lorentz invariance, studying 7^ instead is not equivalent: 
the physics of the dynamics of real plasma excitations is not the same as the physics of 
electric screening. The moral is that defining the Debye mass by the location of the relevant 
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singularity in Fig. ^(6) would be somewhat abstract. 

Another possible method for defining an electric screening length is in terms of carefully 
chosen moments of particular correlation functions. For the purpose of illustration, assume 
we had a correlation function that behaved like 

G(r) = ^ e-"»^^' + ^ e-^''^' . (3.8) 

The first term is the behavior we would have in a TC invariant theory; the second term 
represents the mixing with glueball states due to interactions breaking TC, with e the 
amplitude of that mixing. Now consider defining a correlation length by the ratio of moments 

_ !^r'G{r)dr 

This will yield 

eG = NT)-ix [1 + 0(5//)]. (3.10) 

As long as e is small compared to g"^, this will give a correlation length of order 1/gT 
that, at leading order, matches what we want to call the inverse Debye mass. So one could 
simply define what one means by the Debye mass to be precisely 1/^g- The problem with 
a definition of this form is that it is completely convention dependent; the resulting value 
depends on exactly which correlation function and which moments are used in the definition. 
In contrast, lengths characterizing the exponential long-distance decay of correlations, as in 
our original definition for TC invariant theories, do not generically depend on the details of 
the operators used, other than their symmetries. 

If the Debye mass is so difficult to define at finite chemical potential, what about the 
case of non-relativistic QED? Why do physicists generally have few qualms about discussing 
exponential screening in such plasmas? The reason is a matter of scale and what one 
means by "long-distance." Ref. gives a calculation of the charge-charge correlation in 
such plasmas and finds that it does fall algebraically at very large distance. However, as 
one increases the distance in a variety of physical applications, the correlation first falls 
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exponentially for many e-foldings before finally tapering off in algebraic behavior, and so 
the concept of exponential screening is useful in practice. 

Now consider the case of relativistic gauge theories and our toy example ( |3.8D of a 
correlation G{r). The number of e-foldings over which the first term dominates is only order 
ln(l/£:). Unlike our toy correlation G{r), a real correlation will have additional contributions 
from excited time-reflection odd states, with energies of md+0{g'^T): 

G(r) 0(1) e-['^^^+^(^'^)]'- + ^ E O(^) e-°(^'^)^ + (other junk) . (3.11) 

In the range where the glueball contributions are small, r is still too small to suppress 
these excited states unless gln{l/6) ^ 1. Therefore, an approximate definition of the Debye 
mass, in terms of the intermediate-range fall-off of correlation functions, is not useful beyond 
leading order unless the amount of TC violation is extraordinarily small. In particular, it is 
not useful if the mixing e is simply some power g"^ of the coupling. 



IV. THE 0{g'^T) CORRECTION TO THE DEBYE MASS 

We return now to vectorially-coupled theories at zero chemical potential. If one is in- 
terested only in the 0{g^T) correction to the Debye mass, then it is possible to reduce the 
computation of the Debye mass to a much simpler problem than the extraction of correla- 
tion lengths in a four-dimensional theory with a small periodic dimension and dynamical 
fermions. This simplification will emerge from the successive reduction to equivalent effec- 
tive theories describing longer distance scales, as discussed in section The philosophy 



is similar to that applied by Braaten |[T^ to the expansion of the free energy in powers 
of g. (With more work, it could be extended to handle even higher-order corrections to 
the Debye mass.) The result, to be derived momentarily, expresses the 0{g'^T) part of the 
Debye mass in terms of the perimeter law coefficient of adjoint-representation Wilson loops 
in a three-dimensional pure gauge theory. This relation is particularly nice in that it holds 
regardless of which symmetry channel of the 2+1 dimensional theory has the lowest mass 
time-reflection (or 7?.?) odd excitations. 



First, reduce the problem to an effective three-dimensional theory by integrating out 
modes with non-zero frequency in the periodic direction. If we relabel the z axis as "time," we 
want to know the energy of an Aq, together with its cloud of glue, propagating forward in time 
in 2+1 dimensions. Next, make a further reduction to an effective theory for distances large 
compared to 1/ gT ^ so that the bare Aq can now be considered heavy. The resulting effective 
theory is simply a 2+1 dimensional pure gauge theory (plus irrelevant corrections suppressed 
by powers of g). The non-perturbative contribution of the cloud of glue surrounding the Aq 
is not sensitive to whether the Aq is merely heavy or is infinitely heavy. The propagation 
of the bare Ao can then be replaced by an adjoint-charge Wilson line, exactly analogous 
to the way in which extremely heavy quarks in zero-temperature QCD can be replaced by 
fundamental-charge Wilson lines. The non-perturbative piece of the Debye mass is given by 
the energy of the glue required to screen an infinitely heavy adjoint charge, which can be 
extracted from a numerical lattice calculation of the perimeter-law behavior of large Wilson 
loops. Schematically, 



where mpert is a perturbative contribution to the mass and where Am is extracted from 



in the three-dimensional gauge theory. Note that fermions are completely absent from the 
calculation of Am because they decoupled in the three-dimensional limit. Figure ^ illustrates 
the various stages leading to the relation (|4.1| ), which will be discussed in more detail below. 

At each stage in the reduction, it is necessary to match carefully the effective theories 
onto the original four- dimensional theory, in order to keep track of the perturbative piece 
mpert of ( [4.1[ ) correctly. Fortunately, at the order of interest, the matching is fairly simple 
and straightforward. 

Before we dive into the details of matching, notice that this picture of the Debye mass 
makes the presence of the logarithmic 0{g'^Tln{mQ/ g'^T)) correction found in (|1.5| ), as well 



= ^pert + Am 



(4.1) 




for large loops C, 



(4.2) 
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FIG. 8. Various stages in the reduction relating the 0{g'^T) correction to the Debye mass to the perimeter 
law coefficient of three-dimensional adjoint-representation Wilson loops, (a) The correlation of some pair F 
and G of TZ^ odd operators at large separation. A three-dimensional reduction is performed and only the 
= modes are relevant to what follows. (6) The field is integrated out, generating an propagator 
connecting modified insertions F' ~ dF/dAo and G' ^ dG/dAo. Since Aq is heavy, its propagator may 
be replaced by a straight, adjoint-representation, Wilson line. The long-distance fall-off of the correlation 
measures the energy of glue needed to screen a heavy adjoint charge, (c) If one makes the f direction periodic 
with an arbitrarily large period O, then the minimal screening energy of an adjoint charge is determined by 
the fall-off of an adjoint Wilson line wrapping around the f direction. (This is not to be confused with the 
original Polyakov loop L, which wraps around the small z direction.) (d) The fall-off of the long adjoint line 
with increasing length is the same as the perimeter law fall-off of any large, adjoint Wilson loop. (The loop 
depicted in the figure is meant to be large in x± as well as in f .) 
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as its coefficient, trivial to understand. Consider the self-energy of a static, infinitely heavy 
charge in 2+1 dimensions. If we were discussing QED coupled to such a charge, then the 
electric field surrounding the charge would be £^ = e/27rr, and the energy of that electric 
field would be 

which is logarithmically divergent in both the infrared and the ultraviolet. This same picture 
holds for a non-Abelian gauge theory at distances small compared to the confinement scale, 
since then gluon self-interactions are small. The confinement scale, however, provides an 
infrared cut-off at r~^~(7^T. The mass of the charged particle, which is order mo and not 
actually infinite, provides the ultraviolet cut-off. Finally, since the heavy particle is an 
adjoint charge, should be replaced by Ca g^, where Ca is the quadratic Casimir number 
for the adjoint representation, or for SU(A^). Eq. ([4.3|) then precisely reproduces the 



logarithm of (|1.5|) 



We now turn to fleshing out the details of the split up ( [4.1D and matching of the sequence 
of effective theories. At each stage, an effective theory will describe the same long-distance 
physics as its shorter-distance predecessor, provided its parameters are carefully matched 
to the parameters of its predecessor. One can achieve this matching by computing and 
equating a set of long-distance quantities in both theories. The required matching reflects 
the different treatment of s/ior^-distance physics in the two theories; it doesn't depend in 
detail on the physics at long distances, which in our case is non-perturbative. Hence, if one 
temporarily modifles the theories by introducing a long-distance cut-off, and one uses the 
same cut-off for both theories, then the matching of the infrared cut-off theories will also 
provide the correct matching for the theories when the infrared cut-off is removed. The 
temporary introduction of an infrared cut-off is merely a convenience which allows one to 
compute the matching perturbatively. 
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Step 1: Reduction to three dimensions 



This reduction has been extensively treated in the hterature [pi)t[r4|l4[]^|15| . At the 
order we are interested in, the matching is very simple. Integrating out the non-static 
(Po 7^ 0) components of the fields generates a mass term for Aq through diagrams such as 
Fig. ^. The effective theory is of the form 

83 = ^ J d^x[ + KZ^adj^o)^ + \mlAl + (higer-order) ] . (4.4) 

The "higher-order" term denotes marginal and irrelevant operators in the effective theory 
which are suppressed by explicit powers of g and whose effect can be ignored at the order of 
interest. The dimensionless coupling g is related to the original four-dimensional coupling 

by 

f = g\T)+0{g'), (4.5) 

where evaluating the four-dimensional coupling g &t a. renormalization scale of order T 
eliminates large logarithms from the higher-order corrections to this matching condition. 
The mass mo is just, up to 0{g^T) corrections from two loop diagrams that don't concern 
us, the mass ( |1.4|) we introduced in the introduction: 

< = 1(ca + T. t^) + 0{g'T') . (4.6) 

We have been slightly more general than the SU(A^) case of ( p..4| ); Ca is the quadratic 
Casimir for the adjoint representation and tp is the normalization of each irreducible fermion 
representation F: 

C^6^^ = f^'^dfbcd ^ ^ab ^ (t^^tI) . (4.7) 

We have not yet specified an ultraviolet renormalization scheme for the effective theory 
( [4. 41 ), because it is irrelevant in this step at the order shown above. However we will need a 
specific scheme in the next step, so let us pick one. We will use dimensional regularization 
in d = 3 — 2e dimensions for our effective theory, with the minimal subtraction scheme and 
a renormalization scale fi = T. 
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Step 2: Replacing Aq by a Wilson line 



Now we want to integrate out Aq and move to an effective theory for momenta small 
compared to mo- But our goal is to describe the propagation of an Aq itself, which we 
might probe by the long-distance behavior of some gauge-invariant correlation in the three- 
dimensional theory. For example, consider the three-dimensional analog of the first 0"'"" 
operator listed in Table I, 

(tr(AoFi2)(o,,)tr(AoFi2)(o,o))- (4-8) 

Imagine evaluating this correlation by first doing the path integral over Aq and only later 
doing the path integral over A- The integral over Aq will replace the Aq's appearing in the 
integrand by a propagator of the Aq field in the background of A^ which is the solution to 

^-DfDf + ml) A^„,^(x) = 6%x) . (4.9) 

Now suppose that A has only low-momentum components and is smooth on the scale of 
1/mQ. The solution to ([4.9|) can then be expressed as an expansion in powers of derivatives 
of A and gives 

Am„A{x) = Amoix) Vexp {ij^dx ■ 4adj) + 0{dA) , (4.10) 

where the integration path is a straight line from the origin to x, and where A^o is the free 
scalar propagator which behaves as 

g-mo|x| 

Ar„o(x) ~ — — — - for 1^1 ^oo. (4.11) 

Hence, the net effect of integrating out Aq will be to replace the pair of Aq's in a correlation 
such as ( |4.8| ) by an adjoint-representation path-ordered exponential, times an overall factor: 



^This is true in the approximation that we ignore quartic and higher-order interactions among in 



the effective Lagrangian (^^). As mentioned before, such terms are suppressed by expHcit powers of the 
coupUng and do not contribute to the Debye mass at the order we are interested. To go to higher orders, 
one would have to include these terms and treat them as perturbations. 
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(Ao(x)Ao(O) ■■■)-. e-™il^l ^Pexp(^ dx ■ 4adj) ■ ■ ■) • (4.12) 



This corresponds to Fig. ||(&), and naively mi is just mo. The substitution ( |4.12| ) will be 
valid for separations large compared to the inverse Debye mass. 

But, as sketched in Fig. |^(c), once one introduces the adjoint representation line, one may 
dispense with the details of the operator insertions in the original correlation function by 
considering the f direction (the 2+1 dimensional "time") to be periodic with an arbitrarily 
large period and computing the expectation of a straight adjoint Wilson line wrapping 
around the f direction: 

e-'"^®^trPexp(z^°df -Aadj)^ . (4.13) 

In the Hilbert space interpretation, this corresponds to a trace over all states containing a 
single adjoint-representation external charge. (And includes a sum over all J^^^^ sectors.) 
The exponential fall-off of this expression with the period will be determined by the energy 
of the lightest such state — which is precisely our definition of the Debye mass. 

Finally, as indicated in Fig. §((i), the same coefficient for the fall-off of the correlation 
with contour length may instead be obtained by considering a large, topologically trivial, 
adjoint loop: 



g-mi \C\ 



trVexp^i <j> dx ■ 4adj)^ • (4.14) 



The exponential fall-off of this expression with contour length \C\ will yield the Debye 
screening length through 0{g'^T). The perimeter- law decay of the Wilson loop gives the 
energy of the glue surrounding the heavy adjoint charge, while mi above is the mass of the 
bare charge. 

There is just one complication: the substitution ( |4.1CI| ) is only a good approximation in 



the presence of gauge fields A with small momentum. Our two effective theories — the one 
with Aq and the one where we've replaced it by a Wilson loop — differ in how they treat large 
momentum effects. As usual, this means we need to carefully adjust parameters in order 
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to make the two theories describe the same long-distance physics. In particular, the correct 
choice of mi in ( |4.14| ) is not necessarily mo; it must be determined by matching. This is 
done by perturbatively computing the long-distance fall-off of the Aq correlation functions in 
both theories, which we will only need to do at one-loop order. Dimensional regularization 
will be used as our infrared cut-off. 
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+ 
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FIG. 9. The one-loop self-energy for Aq in the three-dimensional theory. 

In the original three-dimensional effective theory, the fall-off of the Aq propagator is 
determined by the position of the pole. The shift in the pole position due to the one-loop 
self energy can be computed from the diagrams of Fig. ^ and gives a contribution to the 
Debye mass of 



SniQ = CaQ fi 



2,,l+2e 



(27r)3-2^' 



1 (1-26 ) , 2K-p2) 



+ rriQ 



+ q)'^ + mo] 



, fc i^r 2 , 2^ q^ + 2p-q 



q^[{p + qY + ml 



which yields 

5mn = — 



in 



-1 + ln^"^' 



+ 7^-1 



(4.15) 



(4.16) 



We have worked in covariant gauge with gauge parameter ^. 




FIG. 10. The one-loop self-energy for a Wilson line. 



Correspondingly, to obtain the fall-off of an adjoint- charge Wilson line we need to com- 
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pute the one-loop self-energy correction for a static source, as shown in Fig. 0. Locally, we 
can treat a very large Wilson loop as straight. Alternatively, we can confine the system to 
a very large but finite volume and let a straight Wilson line wrap periodically around the 
space, as in Fig. ^{c). In any case, taking the line to be in the z direction. Fig. p!0| gives 

5m, = CAffx dzA,,{z,0) = 0. (4.17) 

^0 

where Aij{z, x±) is the A propagator and the integral vanishes in dimensional regularization. 

Matching requires that mo -|- Sttiq = mi + 6mi. Putting ( [4. 161 ) and ( |4.17| ) together, and 
continuing to choose our renormalization scale in each effective theory to be = T, gives 

mi = mo + 

8n 

The coupling g here still represents the coupling g{T) in the original four-dimensional the- 
ory. The fact that the matching of couplings between effective theories becomes non-trivial 
beyond leading order won't enter into our results at the order of interest. 

Step 3: From the continuum to the lattice 

To measure the perimeter- law behavior of large adjoint Wilson loops numerically, one will 
put the system on a lattice.S The ultraviolet will be regulated by the lattice spacing instead 
of by dimensional regularization, and so we need to modify our matching condition (|4.18|) 
appropriately. We will do this by again matching the one-loop self-energy, now between an 
adjoint line in the continuum and one on the lattice. However, dimensional regularization 
is no longer a good choice for our temporary, common, infrared cut-off in the two theories. 
Instead, we shall consider two, opposite, parallel Wilson lines running in the z direction and 
separated by a large distance R in x. This provides an infrared cut-off because the lines 



"'^^If the lattice theory is defined with hnk variables U in some representation (typically the fundamental 
representation), then the adjoint Wilson loop is given by the path-ordered product of JJadj over the links of 
the loop, where U^^- is tv{T''UT^U^) for each link. 



h In 



7rT2 



+ 7^-1 



(4.18) 
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neutralize each other when viewed from large x±_ = {x, y). The one-loop contribution to the 
energy of these lines is shown in Fig. |1T. 



+ 



+ 



I 



FIG. 11. The one- loop self-energy for two opposing Wilson lines, separated by a distance R. 



In the continuum, the result is 



6'm, = Ca / dz [A^z, 0) - A^z, R)] 

Jo 



(27r)2-2^ ^ 



1 - e'P^-^ 



Stt 



- + ln(7r/i^i?^) + 7£ 
e 



+ 0(6). 



(4.19) 



In lattice perturbation theory, the result is similar but we get lattice propagators instead of 
continuum ones: 



Sm 



lat 



(4.20) 



where a is the lattice spacing; p± is in lattice units and restricted to the Brillouin zone 
\Px\, \Py\ < tt; and 



Aiat(pj 



■ 2 , ■ 2 Py 
sm h sm — 

2 2 



(4.21) 



The a limit of the integral is extracted in the appendix and gives 



2^-1 



ln(8i?V«^) + 27^ 



(4.22) 



To match the theories, we must pick a ^ = fi = T to make the coupling constant definition 
match up, and then require mi^t + 5w^iat = + S'rrii, or 



mi^t = mo + — Ca g^T 
Stt 



8T2 



(4.23) 
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The remaining contribution to the Debye mass is now the quantity extracted from the 
perimeter-law fall-off of large adjoint Wilson loops on the lattice. Since we want specifically 
to extract the coefficient of the 0{g^T) contribution to the Debye mass, we only need the 
leading-order result for the perimeter-law exponent in the limit that g is small. Formally, 
this is 



Am q^T lim lim 

9-^0 |C|->oo 



9'\C\ 



(4.24) 



and |C| is the perimeter of the loop in lattice units. The limit, however, diverges logarith- 
mically as g^O because of the physics behind the logarithm in ( |1.5| ) and ( |4.3| ). To cure 
the problem, we simply need to extract this logarithm explicitly and combine it with the 
perturbative contribution miat: 



I 2rp J , 

md = mo + g T <a + — 

8tt 



In 



rrin 



+ 0{g'T) 



where 



lim lim 



Sg^T'^ J 

j— In <^tr V exp (i j)^dx ■ 4adj) ^ + ^ In (g' 



(4.25) 



(4.26) 



This is our final result for the 0{g'^T) contribution to the Debye mass expressed in terms 
of the perimeter law coefficient for large adjoint-representation Wilson loops in three- 
dimensional pure lattice gauge theory. All that is needed to obtain a numerical result is 
for someone to compute the value of a on the lattice for gauge theories of interest [namely 
SU(3) and SU(2)]. 

We thank Alexei Abrikosov, Lowell Brown, and David Kaplan for useful conversations. 
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APPENDIX 



To do the integral in (|4.20|) , first do the Py integral, which is straightforward and gives 

r, ^ f d^P± A / N /i iv r\ ^ j 1 - cos(2rg) 

Hr) = / 7^ ^Up±) (1 - e^^^^J = — / dq , , (Al) 

J {2ny V ^ 2nJo sin gy^l + sin^ g 
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where q = px/'2- Now split the integral into two pieces: 



lir) 



1 /-Va 



27r Jo 



dq 



1 



singy 1 + sin^ q •? 
When r — s> oo, the first term can be replaced by 

1 



[1 - cos(2rg)] + ^ ^ [1 
27r Jo q 



1 /•'r/2 

2^70 



sing-i/l + sin g 9' 



1 


1 












~ 27 






1 


- In TT + 1 




~ 27 





1 — sin q 



In 2 



The second integral in (^^) is straightforward, and the final result is 

1 



/(r) 



An 



In (8r^) + 27^ 



as r ^ 00. 
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